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The paper is devoted to analysis of different models that describe waves in fluid-filled and 

gas-filled elastic tubes and development of methods of calculation and numerical analysis of 

solutions with solitary waves and kinks for these models. Membrane model and plate model 

are used for tube. Two types of solitary waves are found. One-parametric families are stable 

f^ ' and may be used as shock structures. Null-parametric solitary waves are unstable. The 

(^ ■ process of split of such solitary waves is investigated. It may lead to appearance of solutions 

CN . with kinks. Kink solutions are null-parametric and stable. General theory of reversible 

^ i shocks is used for analysis of numerical solutions. 
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1 Basic equations and formulae 

Solutions of equations of motion for elastic incompressible cylindrical tube with fixed internal 
and external pressure [T] are analyzed here 
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.i_il ' Here ' denotes time differentiation by variable Z, which is Lagrange initial longitudinal 

O ■ spatial coordinate, ' denotes derivation by variable t which is time. Thin membrane model is 

used, unknowns z and r determine surface of the tube in cylindrical system of coordinates, 
z axis of this system coincides with central line of the tube. Parameter p^, is difference 
between internal and external pressure, unknown p is pressure in the material. Parameter p 
is density of material per unit of square so thickness of material which is denoted as h is not 
^_- the parameter of these equations. Parameter h will be essential below in section [3 Here Aj 

are principal stretches and a, are principal Caushi stresses, W is stress function. 
It is assumed that in the case of absence of any external forces 

z = Z, r = R, h = H 

Hence due to condition of incompressibility 

A1A2A3 = 1 

c^i = kWx, 
W = W{Xi,X2A/{Xi,X2)) 
Stress function in this paper corresponds to Gent material 

1 7 1 ^ Xl + Xl + Xl-3 
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In the case of investigastion of solitary waves (section H]) initial data is taken in the form 

z' = z^ + A2'(Z), r = roo + Ar(Z), h = h^o + Ah{Z); Ar,Az',Ah^O,Z^oo 

z{Z) = [ z'iOdC 
Jo 

Condition of equilibrium at infinity is posed 
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There are two equilibrium states at plus infinity and minus infinity in the case of Riemann 
problem (section [5]). Moreover, domains corresponding to other equilibrium states appear 
as a result of solution of this problem. That is why analysis of dispersion relation (section 
[2]) is fulfilled for some arbitrary equilibrium state marked as zero state. 

Notations Wi = W\-^ and so on will be used below. 

2 Analysis of dispersion relation 

In order to analyze results of calculations examination of dispersion relation [2] is neces- 
sary. Dispersion relation may be obtained after substitution oi z = ziexp{i{kZ — uj)t), 
r = riexp{i{kZ — uj)t) into linearized version of equations 01. ip . Note that here physical 
frequency and wavenumber in Lagrange approach is used because all graphs and calculations 
correspond to this approach also. In Euler approach wavenumber equals to k/Xi. It is the 
physically observed wavenumber. Equations are linearized near some uniform state z'^ and 
ro- 
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Let's introduce some notations 



g = Wr,o/{pR), / = aio/(p4'), Ui = ^, U^ = ^/] , 
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li g > f then for sign + before internal square root (plus branch) u/k — )■ Ui when 
k — )■ +CXO. Hence this branch may be treated at infinity as longitudinal branch . For sign — 
before internal square root (minus branch) u/k — )■ Ur when k — )■ +00. Hence this branch 
may be treated at infinity as transversal branch, li g < f then for plus branch u/k — )■ Ur 
and for minus branch u/k — )■ Ui. Examples of graphs of dispersive curves are presented in 
section H] (figl2] and figE]). 

For plus branch + u -^ Uq, c ^ +00 when fc — > 0+ and for minus branch w — )■ 0, c — > f/g, 
Uq is some finite value. 

According to Petrovskii theorem p] equations (11. ip are non-correct if Ui or Ur is imaginary 
value. The first restriction of correctness is natural for membrane model assumption ai > 0, 
the second is also natural for elasticity assumption that stretching of material case resisting 
force. 



The necessary stability conditions: uq ^ R, g > {P*rQ — W2iqY/{pR)- Obvious the 
necessary stabihty condition for uniform state is equihbrium condition 
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But uniform state may be unstable even if equilibrium condition is fulfilled. Such states were 
really found in analysis of solutions of Riemann problem (section [5]). 

3 Solutions of Boussinesq equations and analysis of nu- 
merical methods used in this paper 

This section is devoted mainly to development of numerical methods. Solution of Boussinesq 
equations is treated as test example for development of methods of investigation. 

Boussinesq equations were derived for description of low-nonlinear motions in the paper 
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Here r, ^ and V correspond to t, Z and r — R^o after some stretching. We can also choose 
this stretching so that Ci = 1, because in the case Ci < according to Petrovskiy theorem [3j 
this equation is not evolutionary that means that it is not correct in the sense that Caushy 
problem for such equation is not well-posed. 

Standing solitary wave solution for this equation takes the form [1] 

K(0 = ^sech^ (0 (3.3) 

Note that this equation possesses also travelling solitary wave solutions. Relation between 
standing solitary wave and travelling solitary waves will be revealed below. 

Linearized version of equations (13. 2p was derived in [IJ. Linearization was fulfilled in 
vicinity of solution (13. 3p in order to investigate stability of standing solitary wave. Then 
method of Evans function was used to determine real positive eigenvalue s the square of which 
seems to be the rational number s = 3/16. This fact permits to determine eigenfunction 
analytically 

5(0 = -sech ( - ) + 2sech3 ( - 

Hence solution (13. 3p must be unstable. 

Results of numerical calculations of (13. 2p are described below. The purpose is develop- 
ment of methods of numerical analysis which will be used in section H] for more complicated 
equations. 

Centered three-layer numerical scheme with second order accuracy was used 

Vk\i + Vk% - ^Vk QT' - QV _ n+2 + Vk-2 - 414+1 - 4VT-1 + QVk 

, (^^)g+i + (^^)Li - ^{v'Tk vt' - vr' ^ ^n 
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Here k and n are spatial and time indexes of numerical mesh, A^ and At are spatial and 
time steps, Q = Vr. 

Similar non-dissipative schemes were effectively used in various problems of theory of 
nonlinear dispersive waves such as water waves under elastic sheet waves in composite ma- 
terial and electronic magnetohydrodynamics |1]. The other appropriate scheme is two-layer 



Lax-Wendroff type scheme. For this scheme approximation of spatial derivatives is the same 
as in scheme given above but calculation of time derivatives is the same as in Runge-Kutta 
method for ordinary differential equations. This scheme is also scheme with second-order 
accuracy but it possesses higher-order dissipation. Three-layer symmetrical scheme is more 
preferable in non-dissipative problems and Lax-Wendroff scheme is used for verification of 
results. Additional dissipative terms (viscous elasticity) may be included in both schemes if 
necessary but it is made in different ways. Both schemes are conditionally stable. Condition 
of stability may be obtained by spectral method or by numerical experiment. For A^ —)■ 
it typically takes the form Ar < cA^', c = const. For three-layer method / is usually equals 
to degree of grough co{k) for k -^ oo, here u = uj{k) is dispersion relation of equations un- 
der consideration, or to order of highest spatial derivative. In all calculations in this paper 
relation between time and spatial steps is far from the region of linear numerical instability. 

For non-scalar nonlinear equations slow nonlinear instability may appear if three-layer 
scheme is used even if condition of linear stability is fulfilled. Note that such instability 
was observed for equations with low-order derivatives only. Calculations for fixed-pressure 
equations showed that sometimes nonlinear numerical instability really appears that is why 
when it happened Lax-Wendroff scheme was used. Equations of fluid-filled model considered 
in the last section are similar to gas dynamics equations. Nonlinear instability of three-layer 
scheme for gas dynamics is well-known so only Lax-Wendroff scheme is used. 

Calculations are fulfilled for some bounded domain large enough to avoid essential reflec- 
tions. Rigid boundary conditions (V = const, r = const, z = const) are posed. 

Programs written on Fortran language are used. 

No any graphs are presented in this section because the same results were obtained by 
calculation of full fixed-pressure equations. Similar graphs are presented in the section below 
(figJUandfigJl. 

Firstly initial data for numerical experiment was taken in the form 

V = Vs + eB, Q = esB 

Just as expected for e > on the first stage of the evolution maximal value of r grows 
with time and for e < it decreases. 

In the case e > the grough rapidly increases with time and calculation stops due to 
overflow. This fact may be treated as blowup. 

In the case e < on the final stage of the process standing solitary wave splits into 
pair of two travelling solitary wave moving in opposite directions. No any other waves are 
observed. We can make reverse of the process and treat standing solitary wave as result of 
collision of two solitary waves moving to each other. So standing solitary wave is special res- 
onance solution. Resonances between solitary waves are known for two-dimensional models 
(Kadomtsev-Petviashvili equation) and for models with high-order dispersion (generalized 
Korteweg-de Vries equations and composite material equations) [1]. But in these cases one- 
parametric families of solutions exist. Some of them are stable. The fact that standing 
solitary wave is the result of interaction of two solitary waves suggests to make conclusion 
that such split instability of null-parametric solitary waves will be observed for other models 
typically |5j and that for one-parametric resonance solitary waves such split will be observed 
in some cases. 

Equations 01.11) and 03. 2p belong to class of symmetrical equations for which waves moving 
with some fixed speed are described by systems of symmetrical travelling wave equations. 
Methods for obtaining of symmetrical solitary waves developed in [4j may be applied for 
linearized versions of these equations to obtain symmetrical eigenfunctions. Preliminary 
estimates shows that corresponding shotting procedures in the case of fourth-order equations 
require variation of one parameter if eigenvalue is known from method of Evans function or 
two parameters if eigenvalue is unknown. So calculations seems to be not rather complicated. 



But investigation fulfilled below suggest the method that do not requires calculation of 
eigenf unctions from travelling wave equations for analysis of stability of solitary waves. 
Let initial data is 

V = V,, Q = 

Numerical experiment shows that after some time of standing maximal amplitude groughs 
and then blowup happens. Moreover it was opened that after some time of grough the 
difference B = V{T^:) — Vg is similar to eigenfunction B. Initial data 

V = Vs + eB, Q = eQ 

cause increase and blowup for e > and decrease and split into two solitary waves for e < 0. 

Results of last numerical experiments require some explanations. 

It is well-known from mathematical physics that for linear models with complete count- 
able systems of eigenf unctions (it is typical for problems posed for bounded domains) evolu- 
tion after some time is described mainly by eigenfunction with maximal positive eigenvalue 
because any arbitrary disturbance typically contain some finite component in expansion by 
eigenfunctions. But here we have only one eigenfunction with positive value and problem 
for infinite domain described by dynamical system of continuum type. From formal point 
of view expansion of disturbance for arbitrary initial data typically has zero component of 
this eigenfunction. It seems that here eigenfunction in some sense is function of maximal 
grough during the evolution and all disturbances that are close to this function lead to evo- 
lution that is close to evolution caused by disturbance proportional to eigenfunction. This 
is explanation of observed phenomenons. 

Type of evolution of numerical solution (growth or decrease) in general case depends on 
the difference between numerical solitary wave and solitary wave taken as initial data. Other 
numerical method and other method to obtain initial data are used in section below. Initial 
decrease is observed for e > and blowup is observed for e < 0. 

The other way to get approximation for eigenfunction without analysis of solutions of 
ordinary differential equations is to solve partial equations that are linearized near solitary 
wave. Some initial disturbance is required otherwise zero solution will be obtained. Eigen- 
function of maximal growth (or decrease) will dominate after some time. But in general 
case for some special types of disturbances eigenfunction may not appear. Only oscillations 
will be observed. This is the main difference between linear and nonlinear equations. In 
nonlinear case all waves are related with each other. Hence any arbitrary disturbance leads 
to appearing of all types of waves in nonlinear case. 

Note also that in fact numerical experiments for determination of eigenfunction of max- 
imal growth are fulfilled for some bounded domain. 

4 Calculation of fixed-pressure equations 

The main difference from the previous section is that here standing solitary wave solution 
Tg, Zs is numerical solution of travelling wave equations pQ obtained by the aid of some 
program written on Mathematica language. This program was given by Y.B. Fu. This 
solution depends from some parameters related with conditions at infinity. In the limiting 
case standing solitary wave takes the form of two kinks shifted to large distance one from 
another. According to method of Evans function standing solitary wave solution is unstable 
(exception is the kink case) but no analytic solution for eigenfunction is available. 

Three-layer non-dissipative centered numerical scheme with second order of accuracy is 



given below 
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Instability is observed for long-time calculations. 

Lax-Wendroff type scheme is more preferable here. 
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Calculation with initial data 

r = rg, z = Zg, r = 0, i = (4.4) 

after some standing lead to decrease of maximal amplitude and split of solitary wave. For 
the case of small amplitudes stationary solitary wave splits into two solitary waves travelling 
to opposite directions, figUl So we have just the same type of evolution that described in 
previous section. 

Dispersion curves are presented in figJ2i There is no intersection between dispersive curves 
and line 1 corresponding to speed of travelling solitary waves. Solitary wave of elevation as 
exact solution may exist according to [1] in principal in this case but lower dispersive curve 
corresponds to the case of negative dispersion while Boussinesq equation from the previous 
section corresponds to the case of positif dispersion hence it may describe some virtual waves. 
Line 2 corresponds to speed of linear waves calculated according to results of calculations 
of Boussinesq equation considered in previous section. According to Boussinesq equation 
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Figure 1: Small amplitudes: split 

there is some constant ratio between speed of travelling solitary waves appeared after split 
of standing solitary wave and speed of linear waves. It equals approximately 0.43. This 
line is tangent to minus dispersive curve hence Buossonesq equation describes these solitary 
waves very good. 




Figure 2: Small amplitude: dispersive curves 



The method developed in previous section was used to get solution with increase of 
amplitude. The difference is that here e < 0. 



Calculations with the initial data 

z = Zs + e{z{U) - Zs), r = rs + 6{r{Q -Ts), z = 6z{U), r = er{U) 

give here increase of amplitude are for e < and decrease for e > 0. Here -2;(t*), r(t^:) is 
solution for initial data f l4.4p for some time t* when decrease is already visible. In the case of 
small amplitudes after some time of increase calculations are stopped due to overflow, figl3i 
Blowup may be. So evolution is the same as described in previous section. 
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Figure 3: Small amplitudes: blowup 



Behavior for the case of moderate amplitude differs from behavior in the case of small 
amplitudes, figJH 

In the case of decrease stationary solitary wave splits into two waves also but this is not 
solitary waves. Amplitude of solitary waves decreases due to dispersive effects. Short waves 
are issued. FigJS] shows dispersive curves for this case. Dispersion relation is obtained after 
substitution z = zi exp{i{kZ—oj)t), r = ri exp{i{kZ — u)t) into linearized version of equations 
(11. ip . It was obtained in [2]. Note that here wavenumber in Lagrange approach is used 
because all other graphs correspond to this approach also. In Euler approach wavenumer 
equals to k/Xi. It is the really observed wavenumber. Line 1 corresponds to speed of 
travelling wave measured from numerical experiment. It intersects lower dispersive curve 
for k > 0. According to analysis of well-posed tasks for numerical obtaining of solitary- 
wave solutions [4J there are no travelling solitary waves associated with longitudinal branch 
for equations (11. ip . Travelling solitary waves mentioned in the previous section are only 
approximate solutions here. If initial data will be taken in the form of travelling solitary 
wave the tail of shot waves will appear after some time and amplitude of this wave will 
decrease. Nevertheless special resonance solitary wave exists. Note that by this feature the 
model (ll.ip is similar to model of composite material: classical solitary waves with monotonic 
behavior at infinity do not exist but resonance solitary wave with such property was found 

In the case of increase for moderate amplitudes after some time this increase stops (figED 
and envelop of graph of resulting solution becomes self-similar (fig|7]). It depends oi x/t. This 
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Figure 4: Case of moderate amplitude: initial decrease and spli 
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Figure 5: Case of moderate amplitude: dispersive curves 



solutions contains two solitary-wave like shocks (measured speed of the shock corresponds 
to line 2, there is intersection, chaotic solution may appear after some time) and two kinks 
moving in opposite directions. Analysis of graph for z' shows that short longitudinal waves 
really appear here. Obviously if stationary solitary wave will tend to combination of two 
kinks then speed of kinks in this solution and amplitude of solitary wave-like shock will tend 
to zero. Half of graph these solutions (left or right side from axis of symmetry) coincides with 
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Figure 6: Case of moderate amplitude: initial growth and stop 



the graph of the special solution of Riemann problem, see section below. Typical solution 
of Riemann problem contains two shocks or simple waves and one kink. Here only kink and 
one shock are observed. Hence maximal amplitude r = ro2 may be determined by variation 
of amplitude of initial shock in solution of Riemann problem. 
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Figure 7: Self-similar solution 



For verification of results eigenfunction was calculated also as solution of linearized equa- 
tions, figlHl An interesting fact is that eigenfunction as a result of calculations for a long 



period of time is obtained not for all initial disturbances. For some initial disturbances 
oscillations are observed for a long time. 
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Figure 8: Moderate amplitude: eigenfunction 



Though amplitude is moderate eigenfunction is similar to eigenfunction of Boussinesq 
equations. 

Calculation here are fulfilled for R = 1, p = l,/i = 1, J^ = 30, the same values 
correspond to graphs in sections below; z'^ = 1.1, r^o = 1.69 (small amplitudes), Voo = 1-55 
(moderate amplitudes). 



5 Arbitrary shock split and kink solutions 

In order to verify stability of kink solution Riemann problem of arbitrary shock split was 
solved. Initial data was taken in the step-like form 



roi + {ro2-roi)tanh{{Z - Zo)/L)/2, f = 0, roi 
■4i)tanh((Z-Zo)/L)/2, i = 0, ^ 



^ ~ ^01 + (^02 

Here ro2 is some arbitrary value and zo2 is calculated by formula 
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This formula leads to problem of finding of the root of fifth-order polynomial equation which 
is solved by Mathematika package. Solution with property A = Zq2 



> z'q^ is taken. 



Note 



that there is two such solutions 

bifurcating from zero for A = z, 

evolution to out of correctness calculations. Split into two or three shocks is observed. In the 

case of three shocks one shock is kink mentioned above. But typically it is travelling kink. 



The smaller value is taken below. It means that branch 
02 ~" ^01 is used. The higher value leads after some time of 
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Figure 9: Solutions of Riemann problem close to stationary kink solution. Initial shock was 
located at Z = 0. 



Solution with standing kink is found by variation of parameter ro2. Examples of calculations 
with positif and negative speed of kink are shown below, figjQl 

Analysis of shock problem solutions and dispersion relation for uniform state of right 
side from kink shows that line tangent to minus dispersive curve at the origin of coordinates 
do not intersects it. Dispersion for the origin of coordinates is positive so solitary- wave 
like shock seen on one of the graphs is predicted by theory. Families of solitary waves of 
depression tending to kink solution may exist also. 

Examples of solutions for the case of large enough initial shock are presented here. In 
the case of small amplitude no kink is observed. In the case of moderate amplitude uniform 
state 2 is unstable and calculation of Riemann problem leads to emergency stop. 

Now we can analyze another problem essential from the practical point of view. Let we 
have initial data with two opposite steps: 
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-d < Z < d: 



roi {-d > Z)V{Z > d) 



Large scale approach is assumed here, d is some distance. Steps are smoothed in reality. 
Standing kink solution in this problem corresponds to some critical case. If the speed of left 
kink is negative then such disturbance will increase it's area with time and then the speed 
is positive the area of disturbance will decrease. 

The theory of non-dissipative shocks |1]-[T2] is applicable for this model and the inves- 
tigations made above were fulfilled by the aid of methods developed in this theory. This 
theory is based on analysis of dispersion relation and analysis of dimensions of invariant 
manifolds of solutions of systems of equations of travelling waves numerical modelling of 
invariant manifolds and direct calculations of partial equations. So the same structures as 
in models described in |1]-[I2] exists. Extending non-local shocks with solitary wave struc- 
ture are observed then line U = u/k does not intersect dispersive curve and shocks with 
stochastic structure are observed in the case of intersection. 

The equations with controlled pressure are Hamilton equations. But there is no conserva- 
tion of energy here. Hamiltonian here is not energy. It seems that equations of of travelling 



waves for this model are non-integrable. 

6 Derivation of simplified equations 

Simplified equations are required to describe centered simple waves. In order to derive 
simplified equations let's introduce new notation 

u = z' 

Now equations takes the form 



) — W2{u, r, r') + p^,ru = pRr 



\ y/w^ + r'^ 

—W2{u, r, r') + p^ru = 

Then let's take u from second equation and put it to the first equation and withdraw higher 
order derivatives. Resulting equations takes the form 

^{r",r, ..other.second.order.terms...; r.coefficients.only.) = 0, u = r(r) 

Here $ and F are some functions. 

7 Derivation of equations with bending resistance and 
calculation 

For moderate values of e zones where ui < appear in the solution. This means according to 
Petrovski theorem [3] that in this zones equations (13. 2p are non-correct. The other problem 
is that on the boundary of correct and non-correct zone graphs r{z) and r{Z) take the form 
of broken line. So we have discontinuity of derivative here. 

The model formally permits overlapping and breaking but analysis showed that such 
phenomenons may appear only in the case of compression and hence non- correctness. 

One more problem is that there is no dispersion for short waves in this model {u/k — )■ 
const for k — )■ +oo). It may cause infinite growth of envelop. 

Numerical scheme itself possesses it's own smoothing and correctional properties so cal- 
culation continues some time after non-correct zones appear but later shot-length oscillations 
arise in non-correct zones and calculation stops due to increase of these oscillations. Nu- 
merical solution in non-correct case is usually only approximate solution. Usually there is 
convergence of numerical scheme in correct case if some stability condition holds and there is 
no convergence for any choose of relation between time and spatial steps in non-correct case. 
Hence non-correctness may be revealed numerically also. Petrovskii theorem gives only the 
necessary condition of correctness so then correctness must be verified numerically. 

Obviously non-correctness is related to the fact that membrane model can be used only 
for the stretched case. Plate model which takes into account bending resistance works in 
stressed case also. Bending resistance also obviously prevents the plate from braking that is 
from appearance of discontinuities for r'. Correction of equations is described in section [71 
Note that viscous elasticity model may be used for correction also. 

First of all let's note that the manner of derivation here is differs sufficiently from the 
manner in which equations fl3.2p were derived. Derivation here is based on formulae of linear 



theory of elasticity and converting of formulae of Euler approach into Lagrange form. The 
main goal is condemnation of hypothesis made in previous section that implementation of 
bending resistance may resolve problems appeared in calculations. 

Let's start from well-known Germen-Lagrange formula for stable state of thin unstretched 
and unstressed plate 
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Here pi is pressure acting on the plate, rj is vertical displacement, h is horizontal coordinate, 
h is thickness, E is Yung module, a is Poisson coefficient. Action of internal elastic forces 
caused by bending of the plate may be replaced by action of corresponding pressure. 

For small stretch additional vertical force must be included in left-right side of equations 

m 

dfr = -fxRr""dZ 

o 
Cylindrical form of surface and incompressibility of material a = 1/2, E = 3fi are taken 
into account here. 

Equation for r takes the form 

-cr'"" + [Rcr^^] - ^ + p,rz' = pRr 
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Analysis of dispersion relation shows that now equations are correct independently of whether 
0"! > or (Ti < 0. Calculation with the same parameters as in calculation described in 
previous section was fulfilled, c was treated as a small parameter. No any problems are 
found. 

Solitary wave splits in the case of correction. 

8 Solution of equations for fluid-filled model 

Equations of equations filled by fiuid are given below [2] 

(R(^iti) - Prr' = pRz, [[-cr""]] + [Rcn^) " T^ + P''^' = P^^ 

rz' — r'z + vr' -\ — rv' = 0, Pf{vz' — v'z + vv') + P' = 

Here v is fiuid velocity, pj is fiuid density. Correctional term discussed in previous section is 
included here in double square brackets. 

Dispersion relation for these equations was derived in P]. The main difference from the 
case of fixed-pressure equations is that both branches intersect in point a; = 0, A; = 0. 

It is not clear how to calculate these equations. The one way is replace fiuid gas by 
low-compressible gas. The required modification of hydrodynamic part of equations is given 
below 

[pfz' — f)fZ)r^ + 2pfr{fz' — r'z) + [pjvr'^)' = 

P = P{Pf) 

Equation of conservation of mass is modified here. Equation of state is added. The problem 
is wheather dispersion can prevent breaking of waves in gas if c = 0. 
There are three dispersive curves for these equations. 



The other way is to ehminate pressure and get three equations 
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These equations are solved by Lax-Wendroff type numerical scheme. Due to combined 
time-spatial derivatives this scheme is implicit. System of implicit equations is solved by 
method of iterations. Example of calculation of Riemann problem is given below, figJTO| 
v\t=Q = 0, pf = 1. Kink and solitary wave-like shocks are clearly seen. 

For solitary waves used for shock structures for some given values of physical parameters 
and phase speed U typically one solution exist. Such one-parametric {U is arbitrary pa- 
rameter) families of solitary waves in Hamilton systems are typically stable. Here for some 
special choose of v at infinity these solitary waves may be stationary. There are also other 
null-parametric stationary solitary waves pj] that are similar to waves investigated for the 
case of controlled pressure. They exist only for f = 0. Null-parametric solitary waves in 
Hamilton systems are typically unstable because if energy of solitary wave is slightly de- 
creased there is no way to return it back. Instability of this waves was verified in [llj by 
Evans method. 

The fluid-filled model is the model with conservation of energy. It leads in the case c = 
to integrable system of travelling wave equations and only classical solitary waves and kinks 
are expected. 

Note that the model under consideration due to hydraulic approximation implies slow 
variation of r{Z) hence simplifications may be fulfilled such that no space-time derivatives 
will be. But such simplification obviously leads to decrease in accuracy and appearance of 
non-physical effects. For low-pressure case further simplification may be fulfilled and classical 
or generalized (if bending resistance is included) Korteweg-de Vries equations may be derived. 
Detailed analysis of solutions of generalized Korteweg-de Vries equation is fulfilled in |1] and 
112J 
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Figure 10: Fluid-filled model: Riemann problem 
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